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Inorganic nano-graphene hybrid materials that are strongly coupled via chemical bonding usually present 
superior electrochemical performance. However, how the chemical bond forms and the synergistic catalytic 
mechanism remain fundamental questions. In this study, the chemical bonding of the M0S2 nanolayer 
supported on vacancy mediated graphene and the hydrogen evolution reaction of this nanocatalyst system 
were investigated. An obvious reduction of the metallic state of the M0S2 nanolayer is noticed as electrons 
are transferred to form a strong contact with the reduced graphene support. The missing metallic state 
associated with the unsaturated atoms at the peripheral sites in turn modifies the hydrogen evolution 
activity. The easiest evolution path is from the Mo edge sites, with the presence of the graphene resulting in a 
decrease in the energy barrier from 0. 1 7 to 0. 1 1 eV. Evolution of H2 from the S edge becomes more difficult 
due to an increase in the energy barrier from 0.43 to 0.84 eV. The clarification of the chemical bonding and 
catalytic mechanisms for hydrogen evolution using this strongly coupled MoSa/graphene nanocatalyst 
provide a valuable source of reference and motivation for further investigation for improved hydrogen 
evolution using chemically active nanocoupled systems. 



Sustainable hydrogen production has attracted growing attention, and an advanced catalyst for the electro- 
chemical hydrogen evolution reaction (HER) that results in increased efficiency of this important electro- 
chemical process is urgently required' ^\ Recent work showed that the transition-metal dichalcogenide 
semiconductor M0S2 nanocatalyst, which basically consists of two-dimensional S-Mo-S layers (or 'sandwiches') 
that are stacked to various degrees, is a promising electrocatalyst for the HER"* '°. The catalytic reactivity of the S- 
Mo-S layers is associated with their edges while their basal planes were catalytically inert. Nanosized M0S2 with 
exposed edges should be more active for HER electrocatalysis than materials in bulk forms" ''\ Hence, the 
formation of nanoscale contact between M0S2 and another electronically dissimilar material will open up new 
routes toward novel semiconductor based devices using the advanced catalytic activity of this reduced dimen- 
sional material. 

Graphene, among all kinds of carbon materials, seems to be particularly attractive due to its unique two- 
dimensional honeycomb structure that leads to unusual electronic and mechanical properties and may provide an 
ideal support for semiconductor nanocatalysts'''". Depending on the reaction considered, the role of the gra- 
phene support may be limited to the stabilization of small particles that contain atoms with low coordination 
numbers, or may help to stabilize or activate some of the reactants on the semiconductor itself. Indeed, nanolayers 
of the semiconductor M0S2 supported on reduced graphene oxide sheets have shown advanced electrocatalytic 
and photoelectrochemical properties in the production of H2"'''^. It was found that there is a positive synergetic 
effect between M0S2 and the graphene oxide support, enabling the H2 to be easily evolved from the supported 
M0S2 nanolayers"". Unlike conventional catalysts that have larger sizes, the catalytic performance of semi- 
conductor nanolayers is often controlled by quantum size effects owing to their highly reduced dimensions, 
and is therefore sensitively dependent on the properties of atoms that have lower coordination numbers than they 
would typically have in the bulk. A clear understanding of the nature of the influence of the support on the 
catalytic properties of the supported semiconductor nanolayers using theoretical computation may provide 
insight that leads to new options for nanoscale catalysts. 
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Figure 1 | (a) Two views of the free-standing 3X3X1 M0S2 nanolayer before structural optimization, (b) The top image shows the 7X7X1 graphene 
supercell with M0S2 nanolayer sitting area highlighted as a blue shadow, and on the bottom image the atomic structures of the reduced graphene in the 
shadow area before and after structural relaxation, (c) The optimized structure of the MoS2/graphene nanocontact. The cyan, grey and yellow spheres 
represent Mo, S, and C atoms, respectively. The corner sites are highlighted by dashed rectangles in (a). The missing carbon atoms are highlighted as 
yellow circles and the neighboring carbon atoms with dangling bonds are denoted with numbers in (b). 
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A complete picture of the kinetics of the hydrogen evolution 
reaction (HER) involves, firstly, a primary discharge step (Volmer 
reaction) where bonding of hydrogen to the catalyst + e" + * 
H* occurs (* denotes an active site on the catalysis surface able 
to bind to hydrogen), followed by either an electrochemical 
desorption step (Heyrovsky reaction) or a recombination step 
(Tafel reaction) H* + H* H2. To date, the reaction mechan- 
isms on M0S2 were inconclusive because of the wide range of HER 
Tafel slopes reported'"". The process involving hydrogen bonding 
to the catalyst and the recombination reaction is more inherently 
correlated with the catalyst itself and is an essential step of the 
HER. Knowledge of the H* bonding structure and the reactivity of 
the recombination process on the nanolayered M0S2 catalyst sup- 
ported on graphene sheet is therefore of fundamental importance 
for understanding the reaction that takes place during HER. The 
two types of exposed edges which both have predominantly unsat- 
urated atoms, may provide different bonding geometries with 
hydrogen, resulting in a dual functionality for the hydrogen evolu- 
tion reaction. Besides edges, the kink sites that exist as joints 
between the Mo and S edges may also have some effect on the 
overall catalytic activity. The presence of corners and kink sites 
was revealed by experimental studies, but the role of kink sites has 
not previously been considered^". The mechanism and reaction 
pathways of the HER with M0S2 catalysts without support also 
remained inconclusive. Detailed information about the peripheral 
structures of M0S2 of both free-standing nanolayers and when 
supported on a graphene sheet is thus essential in order to under- 
stand the nature of their catalytic activities. 

In this paper, we analyze from a computational point of view 
the chemical and electronic coupling mechanism between the 
M0S2 nanolayer and vacancy mediated graphene support. In addi- 
tion, the nature of the active sites of the MoS2/graphene nano- 
contact involved in the essential step of the hydrogen evolution 
reaction (HER), that is, hydrogen bonding and recombination, 
was studied. The functionality of the M0S2 nanolayer when a 
graphene base is present turns out to be different for the two types 
of edges, wherein the Mo edge becomes more active in facilitating 
the hydrogen evolution reaction with a lower hydrogen binding 
energy jA£i,ond| and hydrogen recombination barrier compared 
with the free-standing M0S2 nanolayer, but the catalytic effect 
of support on the S edge is the opposite. The calculated results 
provide a motivation for further investigations on improved 
hydrogen evolution using supported M0S2 and provide a valuable 
source of reference of the use of the chemically active nano- 
coupled systems as well. 



Results 

M0S2 consists of layers of S-Mo-S sandwiches. The stoichiometry 
and coordination numbers of the edge atoms can differ from the 
bulk. The most active sites are Mo edges passivated by S atoms that 
have low coordination numbers, providing a sulfide rich envir- 
onment". Although the S-saturation of the Mo edge observed in 
experiments is mostly S-rich, by altering the external environment, 
such as exposure to atomic hydrogen, the S-saturated Mo edge can be 
reduced to pure Mo edge again with missing sulphur atoms'^'. The 
binding energies of sulphur atoms to the Mo edges of MoS2/graphene 
nanocontact and free-standing M0S2 nanolayer were compared, and 
as shown in Supplementary Table 1, the binding energies of sulphur 
atoms at the Mo edge decrease almost 30% when the reduced gra- 
phene support is present. Therefore the defective graphene support 
stabilises a M0S2 nanolayer which is deficient in S. Therefore, in this 
manuscript we focus on the case where a M0S2 nanolayer with pris- 
tine edges is allowed to relax to its minimum energy structure since 
this is what occurs under standard reaction conditions and is there- 
fore an environment that often occurs experimentally. Figure la dis- 
plays the two kinds of bare edge of the M0S2 nanolayer, the sulphur- 
terminated (lOlO) edge and the molybdenum-terminated (lOlO) 
edge, as well as the active S- and Mo-corner sites which are high- 
lighted in the dashed boxes. This is motivated by the frequent obser- 
vation of steps and corners in nanosized M0S2 using high-resolution 
TEM telescope"'". To produce a vacancy mediated graphene sheet, 
two carbon atoms were removed from the graphene sheet, directly 
beneath the M0S2 nanolayer, and the structure was optimised. This 
leads to the presence of a number of dangling bonds on the neigh- 
boring carbon sites, as denoted with numbers in Figure lb, which can 
efficiently chemically interact with the M0S2 nanolayer. Figure Ic 
displays the optimized structure of M0S2 on the vacancy mediated 
graphene surface showing the covalent bond that is formed between 
them, i.e. between the dangling bond of the graphene surface (C2 and 
C4 carbon atoms) and the edge sulphur atom of the M0S2 nanolayer. 
These S-C bonding distances are —(1.76-1.78) A, consistent with 
typical bond lengths in most sulphur compounds. The real applica- 
tions of nanocatalysts supported on graphene are limited as weak 
bonding is usually found, involving the 71 orbital perpendicular to the 
graphene plane^"'. On one hand, the underlying graphene can hardly 
modify the reactivity of supported nanoparticles due to this weak 
interaction. On the other hand, the anchored nanoparticles are 
highly mobile on graphene due to the weak interaction and in some 
cases they easily agglomerate into large particles. A common practice 
in heterogeneous catalysis, controlling the catalytic performance via 
tuning the interaction between the reactive nanoparticles and the 
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underlying support, is to introduce sp^ dangling bonds on the 
graphene surface as a result of the formation of carbon vacancies, 
resulting in vacancy-mediated graphene"^'-^*. Most importantly, 
strongly coupled inorganic nano-graphene hybrid materials stabi- 
lized via this sort of chemical bonding usually present superior elec- 
trochemical performance than the traditional weakly contacted 
counterparts. 

We found that the optimized structure of the M0S2 is quite similar 
whether or not the support is present, and therefore its optimised 
geometric shape in isolated form is not displayed. The optimized 
structure in the presence of the support is shown in Figure lb. The 
M0S2 nanolayer undergoes considerable deformation from the bulk 
structure on structural optimization. As is clear from Figure lb, the 
atomic displacements of the inner S atoms on the (lOlO) Mo edge 
sites during the geometry optimization are significant with extreme 
outward relaxation to the position even beyond the outmost metallic 
Mo neighbors when projected along the z axis. The interacting bonds 
of these displaced S atoms to the inner metallic neighbors is either 
elongated or broken. In contrast, the atomic displacements on (lOlO) 
S edge are small. Hence, the optimized outcomes for each edge. Mo- 
or S-, turn out to be similar to each other with sulfide atoms in the 
outmost exposed positions, which make them the most likely to be 
involved in reaction with hydrogen atoms when they are exposed to 
them. After optimization, the Mo-corner is characterized by two 
more Mo-Mo bonds and results in a four- coordinated Mo atom, as 
shown in Figure lb. Two S atoms appear in the optimized config- 
uration of the S-corner, which are bound to each other and form a 
disulfide linkage triangular unit, which has been reported as a stable 
M0S2 edge configuration in sulfide rich environments based on ther- 
modynamic calculation and TEM observation, i.e.'*'. 

Characterization of the electronic bonding between the M0S2 
nanolayer and graphene support was carried out. As shown in 
Figure 2a, there are several energy levels from —8.0 eV until the 
Fermi level aligned at the same positions in the projected density 
of states (PDOS) of MoS2/graphene, and significant overlap (i.e., 
orbital hybridization) was found between 3p orbitals of S and 2p 
orbitals of C in the wave functions, also indicating a strong covalent 
bond is formed between the nanolayer and the support. The char- 
acter of the charge redistribution between the M0S2 nanolayer and 
graphene sheet can be also evaluated as the difference between the 
electron charges of the M0S2 nanolayer and graphene sheet alone and 
the one formed when they are in contact, and is shown in the inset of 
Figure 2a. The red colour denotes gain of electrons and the blue 
colour denotes electron loss. The contact-induced major charge 
redistribution is not equally distributed among all atoms but mainly 
occurs in the contact region between the M0S2 nanolayer and gra- 
phene support, originating from the vacancy- mediated hybridization 
between the C 2p and S 3p orbitals. 

Next, we compare the total and projected density of states for two 
cases: the free-standing M0S2 nanolayer and the one in nanocontact 
with the graphene support, shown in Figure 2b. From these results, 
one can determine the trend of the charge flow in the M0S2 layer as a 
whole. In both cases, the periphery-induced relaxation breaks the 
symmetry of semiconducting bulk M0S2, leading to the presence of 
metallic states at the Fermi energy level. Note that the states of the 
attached M0S2 nanolayer are similar to that of the free-standing one, 
except for two distinctive changes as pointed out by the arrows in the 
top part of Figure 2b: one is the obvious diminishing of the metallic 
electronic states around the Fermi energy level which are dominantly 
derived from the unsaturated S atoms at the periphery sites; the other 
change is highlighted by the arrow around —8.0 ~ —6.0 eV indi- 
cating occupied states are present when the M0S2 is bound to the 
graphene substrate. Based on the analysis of the projected density of 
states to the atomic sites as shown in Figure 2a, the extra state in this 
lower energy level is associated with the bonding characteristic of the 
C atoms in the support in the attached case. 
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Figure 2 | (a) Orbital-decomposed partial density of states for each atom 
in nanocontacting model. The isosurface of the electron charge density 
difference with the isovalue of 5.0 X 10"' I el A" 'is plotted in the inset with 
the red region denoting the electron gain and the blue region the electron 
loss, (b) Total and projected density of states (TDOS/PDOS) of M0S2 
nanolayer attached on graphene sheet and as the free-standing one as well. 
The distinctive changes are pointed out by the arrows in the top part. 

It has been unequivocally demonstrated that the rate of H2 forma- 
tion on M0S2 catalysts directly correlates with the mean M0S2 par- 
ticle diameter"", suggesting that the active sites for H bonding are 
located at the periphery around the M0S2 particle, which would 
correspond to atoms in corner or edge positions in nanolayers above 
the support interface. To find out whether the decreased metallic 
nature of the states associated with the active periphery sites of the 
M0S2 nanolayer in the supported case affects their catalytic activity in 
hydrogen evolution, we first probed the chemical bonding of hydro- 
gen atoms in the free-standing nanoparticle using the superceU 
model. In the peripheral sites (edges and corners) the cationic Mo 
atoms and anionic S atoms usually have low coordination numbers. 
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Both of these are able to coordinate H, and therefore, they are active (g) 
for initial bonding with H atoms. Our calculations identified the S 
and Mo atoms in the supercell which have low coordination numbers 
as potential active sites on which H atoms can be bonded and recom- 
bine (which were located at the periphery of the M0S2 layer). Binding 
energies of single or two H atoms on a free-standing or supported ^ 
M0S2 nanolayer were calculated as ^ 

A£bond= (-EnH + MoS2|/Graphene] — -BMoS2[/Graphene] — n£H2/2)/n (n=l, 2)(1) -2- 



with £nH-MoS2[/Graphene] being the total energy of the hydrogen 
bonded to the complex with or without the graphene support, 
£MoS2[/Graphene] being the total energy of the corresponding catalyst 
model, and £h2. the energy of the H2 molecule placed inalOX lOX 
10 A^^ cubic box. 

A general principle of catalysis, the Sabatier Principle, is that 
optimal catalytic activity will be achieved on a catalytic surface which 
has intermediate binding energies (or free energies of adsorption) 
with the reactive intermediates^^. It is unlikely that the intermediate 
will be activated by a weak bonding with the catalytic surface, and it 
will be easily contaminated if there is a strong bonding with the 
intermediate. So a compromise between these two extremes with 
an intermediate binding energy is key to an efficient catalytic reac- 
tion. In the particular case of hydrogen evolution, past studies have 
demonstrated that this can be quantified by considering the free 
energy of hydrogen bonding AGbond. and this has been found to serve 
as a reasonable descriptor of hydrogen evolution activity for a wide 
variety of metals and alloys"^' The optimum hydrogen binding free 
energy AGbond should be around 0. The phonon contribution to the 
free energy is typically quite small, <5%, which is negligible in the 
supercell model used. Ignoring entropy terms we can set |AGbond| 
approximately equal to enthalpy |A£bond| and we therefore use cal- 
culated value |A£bond! to judge the catalytic activity for a variety of 
peripheral sites of M0S2 based catalysts, and we will assume in the 
following that the closer |A£bond| is to zero, the better the catalytic 
activity of specified periphery sites. The free energy of hydrogen 
bonded states AGbond was also estimated using the approximated 
method proposed by Norskov by simplifying the entropy energy 
difference to — 1/25^^ , where is the entropy of H2 in the gas phase 
at standard condition, for isolated M0S2 nanolayer and MoS2@gra- 
phene as well (Supplementary Figure 1)'''. However, these differences 
of the hydrogen binding energies AGbond at different edges hardly 
change the predicted relative activities for them by using A£bond- 

In Figure 3a, we have plotted the calculated values of A£bond when 
hydrogen is bound to the peripheral atoms of the isolated M0S2 
nanolayer and M0S2 in contact with the graphene substrate and show 
structural configurations for the latter (H bonded isolated M0S2 
layers show similar structures to the attached one so are not dis- 
played). We found that for the free-standing M0S2 nanolayer, H 
can bind to the S- or Mo-edge and Mo-corner, resulting in negative 
A£bond> as shown in Figure 3a. The H atoms bound to the S-atom 
exposed Mo-edge systems yield very stable structures in which the 
two H atoms are bound to different S atoms at different levels of the 
M0S2 distorted sandwich structure. The S-edge less strongly binds 
the H atoms with a jA£bond| which is closer to zero, suggesting for 
isolated M0S2 particles, the S edge could play a major role on its 
catalytic performance for hydrogen evolution. A stable complex of 
H atoms bound on the Mo corner is formed in which the two H 
atoms are bound to the same corner Mo atom. The main difference 
between the isolated and attached M0S2 nanolayers is that the atomic 
sites at the Mo-edge have a positive A£bond in the graphene attached 
M0S2, with a much lower magnitude than that obtained for the 
isolated M0S2 nanolayer. For free-standing M0S2, the S-edge, with 
its low H binding energy, is the predominant factor which would lead 
to a high activity for the HER of the whole layer. The reactivity of the 
graphene supported M0S2 nanolayer might be a better catalyst for the 
HER because the Mo-edge can contribute to its catalytic activity with 
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Figure 3 | (a) The calculated AEbond (H and 2H) at the peripheral sites of 
MoS2/graphene nanocontact and free-standing M0S2 nanolayer, and (b) 
the optimized structures of the reactants 2H. The H atoms are red. 

a much reduced | A£bond| ■ Although in our study the support does not 
play any direct role in bonding or activating the H atoms, except for 
the role of stabilizing and modifying the M0S2 nanolayer, the S atoms 
on the Mo-edge which have low coordination numbers will be less 
trapping of the H atoms. This could be the result of the fact that some 
metallic states associated with the unsaturated peripheral sites are 
missing after the M0S2 is bound to the reduced graphene support, 
and therefore the M0S2 will more easily release H in the next step of 
electrochemical hydrogen evolution. 

The effect of inclusion of a van der Waals correction on the con- 
clusions was considered for the two-hydrogen-atom binding energy 
at the Mo and S edges of MoS2/graphene nanocontact and free- 
standing M0S2 nanolayer. As shown in supplementary Table 2, the 
introduction of vdW correction alters the absolute value of calculated 
binding energy less than 0.05 eV at the Mo edge in both cases, yet 
leaves the binding energies almost unchanged at the S edge. Hence 
the order of the relative stability of hydrogen bonded each edge is 
preserved when the vdW correction is taken into account. 

Among the different possible MoS2/graphene structures with two 
hydrogens, we also selected the edge sites that weakly bind hydrogen 
and the corner sites that strongly trap hydrogen and have calculated 
the complete reaction path for H recombination and release after two 
atoms are attached to the substrate. The initial state was selected to be 
those shown in Figure 3b-e, and the final state was selected to be the 
local minimum energy configuration that a hydrogen molecule 
placed near the edge or corner would adopt. Activation energies 
£3^4 were then obtained as the energy difference between the max- 
imum energy transition state and the initial adsorption complex 
along the reaction pathway, and reaction energies A£ were obtained 
as the energy difference between initial and the final state, and these 
results are shown in Table 1. 
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In Figure 4, we show the minimum energy profile along the H 
recombination path at the a) Mo and b) S edge sites for the isolated 
and attached M0S2 nanolayer, relative to the energy of the initial 
structure in each case. The detailed atomic configurations of the 
initial state, a number of intermediate states, and the final state are 
also plotted. For two H atoms on the Mo-edge, the global process is 
exothermic with an activation barrier around 0.17 eV for the free- 
standing M0S2 nanolayer and 0.11 eV for the attached nanolayer. 
During the reaction path the first step is a rotation of one H atom on 
the S site followed by a later migration to a high coordination number 
inner Mo atom of the nanolayer. After overcoming an energy barrier 
in the subsequent step, we found both H atoms are situated above the 
nanolayer and come closer to each other. The distance between them 
decreases to 0.75 A in the final state, indicating that an H-H bond is 
formed. The lower barrier of H recombination from the Mo edge on 
attached M0S2 could be due to the nanocontact- induced diminish- 
ment of metallic states on filling the antibonding 2ti* orbital of 
formed H2. We note that the final energy of the H2 on the substrate 
is lower than the separated substrate and H2, which might suggest 
that it remains bound at low temperatures. 

The first stage of the recombination and release of H2 from the S 
edge might be qualitatively considered to be similar to that on the Mo 
edge but the energetics differ. The energy barrier was determined to 
be 0.43 and 0.84 eV, respectively, as the H atom is released from the S 
of the free-standing and attached M0S2 nanolayers. For the isolated 
M0S2 nanolayer, rearrangement to form H2 results in an even more 
stable structure, indicating that this site will not be useful for the HER 
at moderate temperatures. For the MoS2/graphene, there is only a 
small difference in energy between the initial and final structures as 
the H2 is formed, however the activation energy is significant, mak- 
ing this process difficult to proceed. Therefore, the evolution of H2 
from the S edge becomes more difficult due to the increase in the 
energy barrier. Nevertheless in the presence of the graphene, H2 
evolution is enhanced due to the easier evolution path on the Mo 
edge sites, resulting from a decrease in that energy barrier. The excel- 
lent HER activity of the M0S2/RGO hybrid catalyst reported by Li 
et al."" can be understood as Mo edge sites are the active sites mainly 
involved in boosting the hydrogen evolution whereas the S edge 
remains inactive. 

The activation energies for the recombination and release of 
hydrogen were determined to be 1.40 and 2.10 eV, respectively, for 
highly trapped S and Mo corners in the case of MoS2/graphene 
substrate. The hydrogen atoms at the S corner stay bound to different 
S atomic sites initially. As shown in Figure 5, when they start to react, 
one of the H atoms remains bonded to the S corner site while the 
other migrates to a nearby Mo neighbour. The subsequent reaction 
process involves both H atoms breaking from the bonded sites and 
the H2 molecule is assumed to form, but has slightly higher energy 
than the initial structure. From the Mo-corner site, the two H atoms 
are initially bound to the same corner Mo atom. The meta-stable 



product on the potential energy surface, that is 0.22 eV more stable 
than the physically bonded H2, consists of two H atoms located on 
bridge positions near the Mo corner site, with H-H optimized dis- 
tances around 0.78 A. Along the further reaction pathway, the H-H 
bond length starts to decrease to a final equilibrium value of 0.75 A. 

Discussion 

In order to understand and improve the H2 evolution process on a 
M0S2 nanolayer without reducing its high activity, we have investi- 
gated by means of DFT calculations the nature of the active periphery 
sites involved in the H bonding and later recombination and release 
of hydrogen on a series of models in which a graphene support is in 
effect. We found that the coupling with the reduced graphene sup- 
port causes noticeable changes in the shape of the Mo-edge in M0S2 
nanolayer and also a certain degree of charge transfer. The optimized 
Mo-terminated (1010) edge in the M0S2 nanolayer turns out to be 
similar to the S-terminated (lOlO) edge with inner S atoms relaxing 
outward to positions that most readily with the hydrogen atoms, if in 
contact with them. The charge redistribution that results mainly 
occurs in the contact region between the M0S2 nanolayer and gra- 
phene support, originating from the vacancy-mediated hybridization 
between the C 2p and S 3p orbitals. When M0S2 is attached to 
graphene, the metallic states around the Fermi energy level in asso- 
ciation with the unsaturated S atoms at the periphery sites in M0S2 
also decrease as electron transferring to form contact path with 
reduced graphene support. 

We also checked all of the possible active sites on the MoS2/gra- 
phene nanocontact using superceU methods that have been shown to 
be suitable for this type of system. Considering the favourable peri- 
pheral sites, the hydrogen binding energies were used to evaluate the 
catalytic activity, and the activation energy for the essential hydrogen 
evolution step of the recombination process was determined. The Mo 
edge of M0S2 nanolayer in the presence of graphene support clearly 
bound the H atoms less strongly, as was revealed by the small 
|A£bond|. which is enough to modify the energy profile so that the 
following H recombination and release process becomes relatively 
easier with a decreasing activation energy to 0. 1 1 eV. In contrast, the 
same H evolution on S edge becomes more difficult due to a higher 
energy barrier of 0.84 eV for the MoS2/graphene. The functionality 
of M0S2 nanolayer when a graphene base is present is different 
depending on the types of edges, wherein an improved synergistic 
catalytic effect was obtained on Mo edges characterized by lower 
hydrogen binding energies |A£i,ond| ^nd hydrogen recombination 
barrier, but an opposite catalytic effect of graphene support was 
found on S edges. 

Methods 

Here we consider a superceU model of M0S2 having a 3 X 3 in-plane periodicity and a 
single layer on the z axis of the unit cell, see Figure la (note that the structure of the 
isolated nanoparticle is not optimized in this figure). This M0S2 model system was 
studied as an isolated structure and also when supported on a 7 X 7 X 1 graphene 
basal plane repeated in a superceU geometry and with two carbon atoms removed (see 
the resulting optimized geometry in Figure lb). The stacking pattern is not considered 
here as for an adatom or highly anisotropic adsorbant it wiU be of essential import- 
ance, but its effect is smaU for large adsorbants especiaUy of high symmetry. For a 
better comparison, the calculations of the isolated M0S2 nanolayer were done using 
the same superceU employed for the MoS2/graphene systems. A vacuum region larger 
than 1 5 A between vertically repeated models was used, which is large enough to avoid 
interaction between the periodically repeated layers or between the one with bonded 
hydrogen. The M0S2 nanolayer considered has a typical crystaUine hexagonal shape 
with the most reactive (1010) facet exposed for both the Mo-terminated (1010) edge 
and the S-terminated (1010) edges. Unless otherwise stated, the position of aU atoms 
was always fuUy relaxed. The bonding of the hydrogen atoms to both the isolated and 
graphene-supported M0S2 nanolayers was examined. 

All calculations were based on density functional theory (DFT) as implemented in 
the PWSCF code-**. The spin-polarized generalized gradient approximation (GGA) in 
Perdew-Burke-Ernzerhof (PBE) format was included'". The Kohn-Sham orbitals 
used to obtain the electron density were expanded in a plane wave basis set with a 
kinetic energy cutoff of 40 Ryd, and the effect of the core electrons was taken into 
account by means of the Ultrasoft Pseudopotential (USPP) method"'\ Given the large 



Table 1 | The calculated activation Ea^t cind reaction A£ energy (in 
eV) for the hydrogen recombination process on the edge and 
corner sites of the MoS2/graphene nanocontact and free-stand- 
ing M0S2 nanolayer. Activation energy E^et is the energy differ- 
ence between the maximum energy transition state and the initial 
adsorption complex along the reaction pathway, and reaction 
energy AE is the energy difference between initial and the final 
state 



£„et/AE(eV) 


MoS2/Graphene 


M0S2 


Mo-edge 


0.1 1/-0.75 


0.17/-0.72 


S-edge 


0.84/0.20 


0.43/-2.2 


Mo-corner 


2.10/2.10 




S-corner 


1.40/0.31 
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Figure 4 | Calculated energy profile involved in the recombination of 2H on (a) Mo-edge and (b) S-edge of M0S2 nanolayer as supported on graphene 
sheet (red curve) and isolated (black curve) for comparison. The optimized structures of selective images along the energy path for MoS 2/ graphene 
nanocontact are also plotted in the inset. The curves serve to guide the eye. 



size of the supercells employed, calculations were carried out at the T k-point of the 
Brillouin zone. The validity of only adopting a T k-point to describe our supercell 
model was tested using denser k grids of 2 X 2 X 1 and 4 X 4 X 1 . The difference in the 
calculated adsorption energy of single hydrogen atom on the Mo edge site of M0S2/ 
graphene is smaller than 0.002 eV, which validates the suitability of our current 
computational methodology. In optimizing the atomic structure, a conjugate-gra- 
dient algorithm was employed and the force convergence criterion was set to 0.01 
eV/A. 



After the evaluation of the active bonding sites for H atoms and its probable 
recombination path, the next step was to find the activation energy along that path, so 
that kinetics of the recombination reaction could be predicted. The climbing-image 
nudged elastic band (CTNEB) method, widely used for locating transition states, was 
adopted to fmd the minimum energy path (MEP) between a given initial and fmal 
state of the recombination reaction'". An initial path was constructed and represented 
by a discrete set of images of the system connecting the initial and final states. 
Adjacent images were connected by springs, mimicking an elastic band and the 
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Figure 5 | Calculated energy profile involved in the recombination of 2H on (a) S-corner and (b) Mo-corner of MoS2/Graphene nanocontact. The 

optimized structures of initial, transition, and fmal structures along the paths are also plotted. The curves serve to guide the eye. 
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tangent of the path is estimated on each image. An optimization of the band, mainly 
the minimization of the forces acting on images, brings the band to the MEP, through 
which the highest -energy configuration cUmbs uphill to the saddle point. 
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